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Abstract: 

We compute the supersymmetric QCD corrections to the decay distribution of polarized 
top quarks for the semileptonic decay mode ?(t) ^ bl^Vj. As a byproduct, we reinves- 
tigate the SUSY-QCD corrections to the total decay width r{t — > W^b) and resolve a 
discrepancy between two previous results in the literature. 
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1 Introduction 



The dynamics of top quark production and decay will be studied in detail at the Tevatron 
and LHC hadron colliders. Moreover, a possible future linear e+e^ collider will allow for 
precision studies of top quarks, in particular in the threshold region [jl]]. Precise experi- 
mental data will be matched by accurate theoretical predictions, which are possible since 
non-perturbative effects in top quark decays are cut off by the large decay width F ^ 1.5 
GeV. Such investigations may well yield hints to physics beyond the Standard Model, 
since production and decay of top quarks involve very high energy scales. In particular, 
virtual effects of supersymmetric particles may affect top quark production and its decay 
profile [^. Supersymmetric electroweak and strong [Q, ^ quantum corrections to the 
total top quark decay width r{t — > W'^b) have been calculated already some time ago. 
In this article we extend those calculations by considering the SUSY-QCD corrections to 
the fully differential decay distribution of polarized top quarks for the semileptonic de- 
cay mode. From this distribution we can easily derive as a special case the SUSY-QCD 
corrections to the total decay width and compare our result to two conflicting previous 
calculations [@,^. 

Our letter is organized as follows: In section 2 we discuss the calculation of the SUSY- 
QCD correction to the differential decay distribution for ?(t) ^ bl\ and to the total top 
quark width. In section 3 we perform a numerical analysis of our results in terms of 
sbottom and gluino masses, taking into account mixing in the stop sector. Section 4 
contains our conclusions. 



2 Analytic results 

The virtual supersymmetric corrections to the tW^b vertex to order are determined by 
the following SUSY-QCD interaction Lagrangian (where we suppress colour and spinor 
indices of the (s)quark fields and q = t,b): 

Lgqq = V2gsT"q [PtgaqR-PRgaqt] +h.C., (1) 

where ga are the Majorana gluino fields, T" = X"/2 with the Gell-Mann matrices 
and {qi^qR} are the weak-eigenstate squarks that are associated to the chiral components 
Pl,r <? = 5(1 TY5) q of the quarks. The squark mass eigenstates are related to these weak 
eigenstates through a rotation: 

(^j) = ( '^"^Of ^.')^Rl (2) 

\ qi J V -sin9^ cos0^ J \ q^ J 
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Furthermore, we need the contribution of the squarks to the charged current interaction, 
which is given in the mass basis {^1,^2} by: 



-le 



\/2sin0 



w 



RiR),t]d^hW^ 



+ h.c. 



(3) 



Consider now an initial state consisting of top quarks at rest with polarization P. For the 
semileptonic decay 



t{pt) ^b{pt)+r{pi)+vi{Pv 



(4) 



the renormalized amplitude including the SUSY-QCD corrections can be written in terms 
of four formfactors (we neglect lepton masses and the mixing between generations): 



-le 



VlsmQwJ {Pt - Pb)^ ~ + iTwmw 



u{Pv)fPLv{Pi 



X u{Pb){fPL 



l+Fz. + ^(6Zi + 5zi) 



+ fPRpR + — [PlHl + PrHr] \ u{pt) 
mt J 



(5) 



In (^, Z'l^ = 1 + hZ'l" denotes the renormalization constant for the top (bottom) quark 
field, which we fix by imposing on-shell renormalization conditions. This is equivalent to 
the method used in [|5|], where only one renormalization constant for the {t,b) doublet is 
used. In that case an on-shell condition can only be fulfilled by one field, inducing a finite 
wave-function renormalization for the other. Accordingly, we find 



i(5Zi + 6zi)=6Z^-la(m2) 



(6) 



where SZ^ and t[t{mj) are given explicitly in Eqs. (6)-(10) of ref. [^. The form factors 
in Eq. (U) are defined in complete analogy to the corresponding ones in Eq. (3) of [|]], 
except for a relative factor niw /nit in the definition of Hi r. We find complete agreement 
for all formfactors. They are listed explicitly in Eq. (1 1) of [|5|] for arbitrary squark mixing 
angles and masses. Therefore we do not write them down here but only remark that in 
the limit of vanishing Z?-quark mass and no squark mixing the formfactors Fr and Hi are 
equal to zero. 



The phase space Rj, of the final state of reaction 
energies and two angles: 



may be parametrized by two scaled 



dR^ 



m. 



32(271)' 



rdxidxbdxd cos Q , 



(7) 



where Xb = 2Eh /nit, xi = '^Ei / n^t ■ The four- momenta and the polarization of the top quark 
are explicitly parametrized in the top quark rest frame as follows: 
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PI = £,(1,0,0,1), 

pt = £fc(l,O,psin0/ft,|3cose/i), 

Pv = Pt-Ph-Ph 

P = |P|(O,sin0sinx,sin0cos%,cos0), 



where 



^=Jl-4zb/xl COS0 



lb 



XjXt - 2{xi+xt - 1) + 2zb 

XlXb^ 



(8) 



(9) 



with the scaled mass square of the bottom quark Zb = ml/mj. The differential decay rate 
is given by 



1 1 



(10) 



where the sum is taken over the colour and spins of the final state. The fully differential 
distribution for reaction @ reads at tree level: 



lep 



Xi{l -Xi-Zb) 



dxidxbdxd cosQ ^ {I —Xb+Zb — Q^ + 



(l + |P|cos0) 



where 



c — 
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128(271)4 sin^Ojy^ 



with 



m 



m 



(11) 



(12) 



(13) 



Our result for the SUSY-QCD corrections to the semileptonic decay distribution reads: 



dT 



.SUSY-QCD 
lep 



Xl{l-Xl-Zb) 



dxidxbdxd cosQ {\ —Xb+Zb — Q^ 

X {(1 + |P| COS0) Re /i + |P| sin0 [cos^Re /2 + sinxlm /2]} ,(14) 



with 



Xld-Xl -Zh) 



+ 



1 - (1 -Xb + Zb] 



Xl[l -XI -Zb 
1 -Xl 



Xl{l-Xl-Zb) 



[HR + ^bHL]: 



(15) 



Xb^sinQib 
2(1 -xi-zb] 



[{\-xi)Hr + ^{Hi^ + 2Fr)]. 



(16) 
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The function Im/2 can only be nonzero if is larger than nig + rn^^^^, where mj-'^*^' denotes 
the mass of the light stop. We will not discuss this case in the following. 

The SUSY-QCD correction to the total decay rate r{t W^b) can be easily obtained 
from ( p^ in the following way: The narrow width approximation is applied, i.e., one 
makes the replacement 

^ d{{p,-p,f-ml). (17) 



In particular, this fixes the scaled b-quark energy to xt, = I —^ + Zb- The three remaining 
integrations are easily performed. Finally, one has to divide out the branching ratio for 
the semileptonic decay of the W, which is achieved by replacing in ( [TtI ) by r(W^+ 
bl+Vi) = Gpml, / (eVln) with Of = e^jiA^flm^ sin^Giy). The result is: 



pi =rO + rSUSY-QCD ^ ro[l+2ReFz, + Re5Z[ + Re5Z^ 

(18) 



+ 2^ReF^ + 2^Re//L + 2^Re//^ 

Czo Cjo Czo 



where the Born decay Fq rate is given by 

8^271 



r° = S^[(l-^ + ^^)'-4z.]^^'Go (19) 



and 



Go = (l-^)(l+2^)+z,(z, + ^-2), 
Gi = ^[(l-^)2 + z,(z,-2^-2)], 

G3 = ^G2. (20) 



Our result for the total decay rate disagrees with the corresponding result given in Eq. 
(15) of [|]]. The disagreement appears to be due to an error that occured in deriving G23 
from the standard matrix elements M2,3 given in Eq. (13) of [|^]. The result in can 
be corrected by interchanging G2 ^ G3 (or, equivalently, Hj^ ^ Hr). In an earlier work 
[0], the supersymmetric QCD contributions to the top quark width have been computed 
for the special case of degenerate SUSY masses and = 0. We performed a numerical 
comparison with Figures 2 and 3 of [Q] and find complete agreement when using the same 
input parameters. 
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3 Numerical analysis 



In this section we discuss the impact of the SUSY-QCD corrections on the total top quark 
decay width, on the energy spectra of the charged lepton, and on observables sensitive to 
the top quark polarization. 

We start by considering the relative correction 

pl_pO 

S| = ^^ (21) 

to the total decay rate. As mentioned above, this quantity has been studied in the literature 
before with two different results. Our calculation confirms the earlier result [Q]. The 
effects of the mixing of the chiral components of stop and sbottom have been considered 
only in [|]]. Therefore it seems worthwhile to reconsider the quantity 5|. 

The stop and sbottom mass matrices can be expressed in terms of MSSM parameters as 
follows: 

^2 _ ( M| + m2 + m|(i-a4)cos2|3 m,(Af - ^/cotP) 

' ~ \ m,(Af-A/cot(3) M^ + mj + mlQtS^cosl^ 



= ( Ml + ml-ml{^ + QbS^)cos2^ mfe(Ai -A^tanp) \ 

^ \ mh{Ab-fitm^) M?+m2 + m|2fo4cos2|3 y ' ^ 

where Mg, Mq.M^) are the soft SUSY-breaking parameters for the squark doublet qi 
and the squark singlets tR and bR, respectively. Further, A^ j, are the stop and sbottom soft 
SUSY-breaking trilinear couplings, and ^ is the SUSY-preserving bilinear Higgs coupling. 
The ratio of the two Higgs vacuum expectation values is given by tanP, Qt = 2/3 and 
2^ = — 1/3 denote the electric charges of t and b, and sw = sinGiy. The squared physical 
masses of the stops and sbottoms are the eigenvalues of the above matrices. In order to 
keep the numerical discussion tractable, we make the following simplifying assumptions: 
We neglect mixing in the sbottom sector. This is certainly justified if tan|3 is not too large. 
In any case tan (3 only enters through the mass matrices. If sbottom mixing is neglected, 
the dependence on tan (3 is very weak [|]] and we set tan (3 = 1 for all following results. 
Further, we set Mq = Mfj and neglect the bottom quark mass in the mass matrices. Under 
these assumptions the sbottom mass matrix is diagonal with degenerate mass eigenvalues, 
M-' = diag(m? , m~ ) . Note that using degenerate sbottom masses close to the experimental 
lower mass limit maximizes the impact of the SUSY-QCD corrections. The stop mass 
matrix simplifies under the above assumptions to 

9 / mr+m? nitMiR \ 

= ,. ,.2 2 , (23) 
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100 200 300 400 

m^[GeV] 

Figure 1 : SUS Y-QCD correction 5g as a function of the gluino mass for different sbottom 
masses and no mixing: = 80 GeV (a), 120 GeV (b), 160 GeV (c) and 200 GeV (d). 



with MiR = A; — /J. Maximal mixing (Bf = | and M^r ^ 0) corresponds to Mq = rn^. 
The latter relation will also be assumed for Mir = 0, leading to the following stop mass 
eigenvalues^: 

m,- 2 = sJml + m}±mtMLR. (24) 

Fig. 1 shows 5| for Mir = as a function of the gluino mass for different values of 
m^. The SUSY-QCD corrections are negative and of the order of several permill for 
gluino masses larger than 100 GeV. Even for very small gluino masses the SUSY-QCD 
corrections are at most ~ (—1)%. Our Fig. 1 corresponds exactly to Fig. 2a of [^. 
In particular, we use nit = 174 GeV and CLs{mt) = 0.11. (For the bottom quark mass 
we use Mb = 4.75 GeV.) We find about 30% to 40% smaller SUSY-QCD effects than the 
authors of and can exactly reproduce their curves if we, just for this purpose, substitute 
Hl^Hr. 

The effect of mixing is studied in Figs. 2a,b, where we plot 5| as a function of the mixing 
parameter Mir for different sbottom and gluino masses. For M^r = 200 GeV, nig = 150 

GeV and = 100 GeV (which implies mf^^^ = 74 GeV), the SUSY-QCD corrections 
reduce the total top quark decay width by about 2%. Larger squark and/or gluino masses 
lead to smaller SUSY-QCD corrections. Note that the squark masses we use are compat- 
ible with bounds obtained in a recent ALEPH analysis [^. For the gluino mass, experi- 

'Note that by fixing Qf = j the light stop can be either t\ or f2 depending on the sign of Mlr. 
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mental lower mass limits are typically higher than 200 GeV (see, e.g. |^]), but these 
limits only apply within the minimal supergravity model. 




Figure 2: SUSY-QCD correction 5| as a function of the mixing parameter Mir for = 
100 GeV (a) and = 120 GeV (b). The full curve is for lUg = 150 GeV, the dashed curve 
for m, = 200 GeV 



We now turn to the discussion of the fully differential leptonic decay distribution, Eq. (P^. 
In Fig. 3a we plot the charged lepton energy spectrum JFiep/ dxi and in Fig. 3b the relative 
SUSY-QCD correction 



Slep(-^z) 



lep 

dxi 



-1 



dV^ 

lep 



dxi dxi 



(25) 



Li the narrow width approximation for the W boson we have 



5iep(^/) =Re/i. 



(26) 



We consider here the case of maximal mixing with Mm = 200 GeV and masses = 100 
GeV and rtig = 150 GeV. In this case 6iep(x/) reaches values of —2.7% close to the sharp 
drop of the energy spectrum at xi ~ 0.2. As can be seen in Fig. 3b, the narrow width 
approximation for the W propagator works well in almost the whole kinematic range for 
xi which is allowed within this approximation. 

A sample of highly polarized top quarks (which can be produced at a linear collider with 
polarized beams operating close to the ti production threshold) would allow for additional 
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Figure 3: SUSY-QCD corrections to the charged lepton energy spectrum: (a) shows 
JFiep/Jx/ in GeV in leading order (dashed line) and including the SUSY-QCD correc- 
tions (full line), (b) shows the relative correction 6iep(jc/) in percent. The dashed curve in 
(b) shows 5iep(jc/) using the narrow width approximation for the W propagator. All curves 
are for = 100 GeV, m~g = 150 GeV, and MtR = 200 GeV. 



tests of the top quark decay profile. A well-known characteristic of semileptonic decays 
of polarized top quarks is the factorization of the double differential cross section 



dxid cosQ 



/(x/)(l + |P|cos0), (27) 



which holds true not only at the Bom level, but also to high accuracy including QCD 
radiative corrections This means in particular that the charged lepton is the perfect 
analyser of the top quark spin, i.e. the distribution (F)^^ JF/JcosG has maximal slope |P| 
up to permill QCD corrections. As exhibited by Eq. ([T4|), SUSY-QCD corrections respect 
the factorization ( PTj ) exactly. This means that the normalized distribution 1 /FJF/JcosG 
is not affected by the SUSY-QCD corrections. 

The general decay distribution (|T^ contains a further term for nonzero top quark po- 
larization, which is determined by the function Re /2. This term may be accessed by 
considering the azimuthal asymmetry 



' P > / . / / ! , • (28) 



I dxi 



J3n/2 Jn/2 J dxidX 
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Figure 4: Azimuthal asymmetry 55^(x/) for |P| = 1 and the same choice of mass parameters 
as in Figs. 3a,b. The dashed line shows 65^ (;c/) in the narrow width approximation. 



Note that ^xi^i) ■^^^^ leading order. In the narrow width approximation, 

IPI 



(29) 



Fig. 4 shows ^xi^i) same choice of mass parameters that have been used in 

Figs. 3a,b and for maximal top quark polarization |P| = 1. The asymmetry is negative 
and of the order of a permill. 



4 Conclusions 



The results of our analysis of the SUSY-QCD corrections to the decay ?(t) — ^ blv may be 
summarized as follows: 

1 . The total decay width of the top quark is reduced by a few permill (no mixing) up to 
several percent (maximal mixing in the stop sector, sbottom masses around 100 GeV and 
gluino masses in the range 150 to 200 GeV). A conflict between two previous calculations 
[Q, ^ has been resolved in favour of the earlier work [Q] . 

2. The SUSY-QCD corrections to the energy spectrum of the charged lepton reach values 
of almost —3% for maximal mixing. 

3. Observables that are sensitive to the top quark polarization are hardly affected by the 
SUSY-QCD corrections: The tree level factorization of jriep/((ijc/c^cos0), cf. Eq. 
is respected, and the azimuthal asymmetry (EHl) induced at one-loop is tiny. 
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